Navigating surface reconstruction of spinel oxides for electrochemical water oxidation

Understanding and mastering the structural evolution of water oxidation electrocatalysts lays the foundation to finetune their catalytic activity. Herein, we demonstrate that surface reconstruction of spinel oxides originates from the metal-oxygen covalency polarity in the MT–O–MO backbone. A stronger MO–O covalency relative to MT–O covalency is found beneficial for a more thorough reconstruction towards oxyhydroxides. The structure-reconstruction relationship allows precise prediction of the reconstruction ability of spinel pre-catalysts, based on which the reconstruction degree towards the in situ generated oxyhydroxides can be controlled. The investigations of oxyhydroxides generated from spinel pre-catalysts with the same reconstruction ability provide guidelines to navigate the cation selection in spinel pre-catalysts design. This work reveals the fundamentals for manipulating the surface reconstruction of spinel pre-catalysts for water oxidation.

core-shell structure, in which the oxyhydroxides in the surface and near-surface region are considered as the real catalytic species for the reaction. This paves the way for new design paradigm to manipulate the surface reconstruction of pre-catalyst for advanced OER activity. However, to date studies about surface reconstruction mainly focus on exploring the compositional effect on the activity of the reconstructed oxyhydroxides; 11,12 it remains unclear in which way and to what extent we can control the surface reconstruction of pre-catalyst.
A prerequisite to manipulate surface reconstruction is to identify what triggers the process of surface reconstruction. Some pioneer works recommend cation leaching as the origin of surface reconstruction 8,13,14 . For example, lithium cation is considered unstable in alkaline conditions and therefore can easily leach out to activate surface reconstruction. J. Wang et al. have achieved severe in situ reconstruction in layered LiCoO 2-x Cl x (x = 0, 0.1 and 0.2) due to lithium leaching 14 . However, no such reconstruction was observed in spinel LiMn 2 O 4 , where lithium cations are intensively involved but the spinel framework remains almost unchanged during OER 15,16 . The contradictory phenomena suggest that the leaching of metal cations is not element-dependent but should be material-dependent, which necessitates deeper investigations into the electronic structures to give an explicit explanation. Unfortunately, this part has been less touched and remains elusive. Besides, precise control of the reconstruction ability of pre-catalyst is demanded. Surface reconstruction may not always bring activity enhancement. If the intrinsic activity of pre-catalyst is superior to that of the reconstructed counterpart, strategies should be made to suppress the reconstruction to guarantee a long reaction durability. On the contrary, if the reconstructed oxyhydroxide shows superior activity, the pre-catalyst should be designed with high reconstruction ability to promote surface reconstruction. This calls for the establishment of a structure-reconstruction relationship, based on which the reconstruction ability of pre-catalyst can be forecasted and finetuned. In addition, principles are needed to guide the selection of metal cations in the design of pre-catalyst. The activity of the reconstructed oxyhydroxides is dependent on the type and ratio of metal cations 11 . Therefore, to facilitate in situ generation of highly active oxyhydroxides, investigations are intensively needed to pinpoint the active cation species in reconstructed oxyhydroxides.
Herein, using spinel Li x Co 1-x Co 2 O 4 as model catalysts, we demonstrate that surface reconstruction of spinel oxides originates from the biased metal-oxygen covalency between MO 4 and MO 6 . Combining computational and experimental approaches, a structurereconstruction relationship has been established, based on which the reconstruction ability of spinel pre-catalyst can be effectively forecasted. Further electrochemical analysis and computational calculations deliver some useful clues for the cation selection in spinel precatalysts. This work unveils the structural origin of reconstruction for spinel oxides. It also provides a structure-reconstruction relationship applicable to all spinel family to forecast their reconstruction ability. The demonstrated strategies to master surface reconstruction provides firm theoretical foundations to navigate the design of spinel precatalysts with controllable reconstruction ability for electrochemical water oxidation.

Results and discussions
Theoretical investigation of the structural origin of surface reconstruction on spinel oxides The structural composition of spinel oxides is AB 2 O 4 , in which the cations in A sites and B sites are tetrahedrally and octahedrally coordinated, respectively. Within the structure, the octahedral units connect to one another via an edge-shared form, while the connection between tetrahedral units and octahedral units is corner-shared (Fig. 1a). The oxygen ligands are arranged in a cubic close-packed lattice and each of them is shared by three octahedral cations and one tetrahedral cation to create a M T -O-M O backbone 17 . When catalyzing OER, the spinel frame may collapse and reconstruct to amorphous oxyhydroxides during electrochemical cycling. The evolution from spinel to amorphous oxyhydroxides necessitates the breakage between the tetrahedral and octahedral units, during which highly reactive and transferrable oxygen ligands are believed to facilitate the phase transformation. Besides, because the cations in the reconstructed oxyhydroxides stay in an octahedral coordination environment, a stronger metal-oxygen covalency of the octahedral units than the tetrahedral units may be advantageous for the reconstruction process. In other words, a biased M T -O-M O backbone with oxygen binds more strongly to the octahedral cations may be decisive for more smoothly and thoroughly reconstructed oxyhydroxides.
To verify above hypotheses, a case study is performed on spinel (Co 1−x Li x )Co 2 O 4 (x = 0, 0.25, 0.5, 0.75 and 1). Spinel (Co 1−x Li x )Co 2 O 4 is chosen because the CoO 6 units are the only catalytically active centers after reconstruction, which allows us to explicitly analyze the degree of surface reconstruction. The density functional theory (DFT) calculations are carried out to explore the electronic structure evolution with the change of spinel composition. As a prerequisite for structural collapse, the reactivity of oxygen ligands is first examined. With the increase of lithium in spinel (Co 1−x Li x )Co 2 O 4 , the oxygen p-band center gradually increases and moves closer to the Fermi level ( Fig. 1b), indicating the oxygen ligands are becoming more active. Meanwhile, the calculated oxygen vacancy formation energy (E vac-O ) of CoCo 2 O 4 is as high as 5.20 eV, suggesting severe difficulty of generating an oxygen vacancy in the structure. With the tetrahedral cobalt cations being replaced by lithium cations, the energy required to generate oxygen vacancy steadily decreases (Fig. 1c). The value of E vac-O decreases to 0.80 eV in LiCo 2 O 4 , implying high mobility of oxygen ligands to escape from the spinel lattice. The trends of oxygen p-band center and oxygen vacancy formation energy demonstrate an activation effect to the oxygen ligands with the increase of lithium in spinel ( The M-O bond strength of the M T -O-M O backbone is subsequently evaluated by calculating the metal-oxygen covalency in the spinel structure. The metal-oxygen covalency is defined as the difference in band center between oxygen p-orbitals and metal d-orbitals, in which a higher value represents a weaker covalency 17 . As shown in Fig. 1d, the metal-oxygen covalency of tetrahedral units and octahedral units are quite close to each other in spinel CoCo 2 O 4 and (Co 0.75 Li 0.25 ) Co 2 O 4 . In this case, oxygen binds equally to the tetrahedral and octahedral cations, which fabricates a nonpolar M T -O-M O backbone. The spinel frame is therefore stable and should be difficult to collapse for reconstruction. However, with more lithium cations being introduced into the spinel lattice, the metal-oxygen covalency of tetrahedral and octahedral units begins to get weakened and strengthened, respectively. In spinel (Co 0. 5  , the calculated covalency gap is as high as 11.03 eV, implying the spinel frame is highly unstable and the reconstruction to octahedrally coordinated oxyhydroxides is considerably feasible. To further investigate the trend of surface reconstruction ability in spinel (Co 1−x Li x )Co 2 O 4 , the surface transformation energy from spinel to layered oxide is calculated. The surface transformation energy is defined as the energy of (111) spinel surface relative to that of the layered CoOO species (more details are summarized in Methods). As shown in Fig. 1e, the surface transformation energy of spinel CoCo 2 O 4 and (Co 0.75 Li 0.25 )Co 2 O 4 is far above 0 eV, suggesting the transformation from spinel to layered species is difficult. With more lithium cations substituting cobalt cations, the energy (per surface area) needed to transform spinel to layered structure gradually decreases. This indicates the difficulty for surface reconstruction is reduced. In spinel (Co 0. 25 Li 0.75 )Co 2 O 4 and LiCo 2 O 4 , the surface transformation energy is quite close to and even smaller than 0 eV, hinting very easy-going surface reconstruction. Therefore, we can rationally deduce that the increase of lithium cations in spinel (Co 1−x Li x )Co 2 O 4 will result in enhanced surface reconstruction ability.

Synthesis and structural characterizations of spinel Li x Co 3−x O 4
To verify the reconstruction trend predicted by DFT calculations, spinel Li x Co 3−x O 4 (x = 0, 0.25, 0.5, 0.75, and 1) were prepared by a nitrate decomposition method followed by heat treatment (more details are summarized in Methods). The inductively coupled plasma optical emission spectrometry (ICP-OES) results (Supplementary Table 1) confirm that all the as-prepared samples remain stoichiometric Li x Co 3−x O 4 . The Li 1 s XPS spectra of the pristine samples and the derived lithium molar ratio are summarized in Supplementary Figure 1. The results evidence that, although slightly leached, the surface lithium concentration still obeys the same trend as the total lithium concentration. The crystal structure of the as-prepared spinel Li x Co 3−x O 4 was confirmed by X-ray diffraction (XRD) analysis. As shown in Fig. 2a, the diffraction peaks of the as-prepared samples match well with that of the standard Co 3 O 4 (PDF#00-042-1467), indicative of typical cubic spinel structure. The peak at 19°gets intensified with the increase of lithium in spinel Li x Co 3−x O 4 , suggesting the gradual decrease of tetrahedral cobalt cations and that the introduced lithium cations mainly reside in the tetrahedral sites 18 . The X-ray photoelectron spectroscopy (XPS) was performed to further evaluate the surface Co 2+ /Co 3+ ratio of the samples. The Co 2p spectrum exhibits two major peaks at 795 eV and 780 eV, ascribed to the typical Co 2p 1/2 and Co 2p 3/2 orbitals, respectively. An approximate energy variation of 15 eV is observed between the two major peaks, consistent with that reported in previous studies 19,20 . The deconvolution of the two major peaks shows two characteristic peaks, which can be ascribed to Co 2+ and Co 3+ , respectively. As displayed in Fig. 2b, with more lithium cations substituting cobalt cations, the relative intensity of Co 2+ characteristic peak gradually decreases, indicating the decrease of Co 2+ cations and an increase of cobalt oxidation state. The quantitative atomic ratio of Co 2+ /Co 3+ could be estimated by comparing the area of the characteristic peaks 19 . The results, which are given in Supplementary Table 2, show an obvious decrease of the surface Co 2+ /Co 3+ ratio, further hinting the occupation of lithium cations in the tetrahedral sites.
The Fourier Transform Infrared (FTIR) spectroscopy was employed to probe the local chemical bond change of the as-prepared samples (Fig. 2c). The bands in the frequency range from 500 cm −1 to 700 cm −1 are attributed to the vibrations of the CoO 6 octahedra 21 . Specifically, the one at around 570 cm −1 is attributed to the asymmetric stretching modes of CoO 6 and the one at around 670 cm −1 to the bending modes of Co-O-Co 21 . The band at around 448 cm −1 corresponds to the stretching vibrations of LiO 4 tetrahedra 22 . This characteristic band gets intensified with the increase of Li in spinel Li x Co 3−x O 4 , which indicates, once again, the residence of lithium cations in the tetrahedral sites. The cation occupancy was also confirmed by the extended X-ray absorption fine structure (EXAFS)   technique. Figure 2d displays the Fourier-transformed k 3 -weighted Co K-edge EXAFS spectra of the as-prepared samples. The peaks at 1.5 Å and 2.5 Å correspond to Co-O scattering and Co Oct -Co Oct scattering, respectively. The peak at 3.0 Å corresponds to the Co Tet -Co Oct scattering. With more lithium substituting cobalt, the intensity of Co Tet -Co Oct scattering peak steadily decreases, suggesting the gradual absence of tetrahedral cobalt in spinel Li x Co 3−x O 4 and therefore implying the occupation of lithium cations in tetrahedral sites. It can be rationally concluded that the as-prepared spinel Li x Co 3−x O 4 .
Electrochemical analysis of the pristine and fully cycled spinel After confirming the crystal structure and cation occupancy of the asprepared samples, the electrochemical analysis was conducted. The electrochemical measurements were carried out within a potential window of~1.0-1.62 V (vs. RHE) in 1 M KOH, with the catalytic activities being normalized by the surface area of oxides (more details are summarized in Methods). The surface area of oxides was obtained by Brunauer-Emmett-Teller (BET) technique, which has been proved to be suitable for measuring the surface area of metal oxides in powder form 23 . The current efficiency of spinel LiCo 2 O 4 is 95.4 ± 2.5% (Supplementary Figure 3), confirming that the measured current mainly originates from OER. Figure 3a summarizes the iR-compensated activity of the pristine samples. The non-iR compensated results are shown in Supplementary Figure 4. The activity rises with more lithium cations substituting cobalt cations, probably because of the strengthened Co Oct -O covalency 24 . As the electrochemical cycling continues, the surface and near-surface region of the spinel structure may collapse and evolve into amorphous oxyhydroxides (Fig. 3b). The structural information and OER activity of the fully cycled spinel Li x Co 3−x O 4 were analyzed after performing the electrochemical cycling for 500 times. The XRD and FTIR results of the cycled samples are shown in Supplementary Figure 5, which evidences that the bulk of the samples remains spinel crystal and implies that the reconstruction process mainly occurs on the surface region. To probe the elemental composition of the near surface region after reconstruction, the Co 2p and Li 1 s XPS spectra of the cycled samples are recorded (Supplementary Figure 6). As illustrated, no Li 1 s spectra can be recorded, suggesting the absence of lithium in the near surface region. Besides, of all the samples, the positions of Co 2p characteristic peaks are quite close to those of CoOOH, which further evidences the formation of surface cobalt oxyhydroxide. Therefore, we believe the active sites after reconstruction should be the CoO 6 units. The O 1 s XPS spectra were recorded to investigate the components of various oxygen species in the surface region. As exhibited in Fig. 3c, the relative ratio of lattice oxygen and surface oxyhydroxide becomes much smaller and higher, respectively, with cobalt being gradually substituted by lithium. The numerical ratios can be extracted by calculating the area the characteristic curve covers, which are given in Fig. 3d. The results evidently demonstrate the steady growth of surface oxyhydroxide with more lithium being present inside the spinel lattice, suggesting higher surface reconstruction degree with the increase of lithium cations. The comparison of O 1 s XPS spectra between the pristine and cycled samples are shown in Supplementary Figure 7. As demonstrated, compared to the pristine samples, all the cycled ones possess increased ratio of surface oxyhydroxide and decreased ratio of lattice oxygen, suggesting the destruction of spinel structure and the in situ formation of oxyhydroxides. Besides, more lithium incorporation is found to promote the generation of more surface oxyhydroxide, which further evidences the intensified surface reconstruction in high lithium-containing Li x Co The change in OER activity gives a more precise evaluation of the surface reconstruction degree. As above-mentioned, all the prepared samples share the same catalytically active centers after reconstruction, i.e., the reconstructed CoO 6 units. This lays the foundation that the activity change between reconstructed catalyst and the pristine one can be convictive for assessing the surface reconstruction degree. The cyclic voltammetry (CV) curves of all the pristine and cycled samples are shown in Supplementary Fig. 8. It is obvious that the introduction of lithium is beneficial for a more apparent catalytic enhancement, as Co 3 O 4 exhibits almost unchanged CV curves while LiCo 2 O 4 shows dramatic rise after fully cycling (Fig. 3g). It's worth to note that, based on the chronoamperometry results of spinel Co 3 O 4 and LiCo 2 O 4 ( Supplementary Fig. 9), the incorporation of lithium in   h The experimentally measured current density (@ 1.62 V vs. RHE) of the pristine and cycled spinel Li x Co 3−x O 4 . i The plot integrating both the experimentally measured and theoretically calculated reconstruction degree. The current density (@ 1.62 V vs. RHE) change (Δ I ) between cycled and pristine spinel Li x Co 3−x O 4 is used to represent the experimentally measured reconstruction degree, while the DFT calculated metal-oxygen covalency polarity (the metal-oxygen covalency difference between MO 4 and MO 6 , Δ C ) is used to represent the theoretically predicted reconstruction degree. The insert shows the linear fitting of the two parameters.
spinel structure does not influence the OER durability, as both catalysts exhibit rather stable catalytic performance within 60 operating hours. Figure 3h provides a direct comparison of the current density (@ 1.62 V vs. RHE) between the pristine and the cycled one. Similar to that illustrated in Fig. 1d Fig. 10). As exhibited in Supplementary Fig. 10a, the increase of ECSA gets intensified with the increase of lithium in pre-catalyst. This reveals the increased active sites, and therefore demonstrates the higher amorphization degree, in high lithium-containing samples. However, the increase of surface area does not alter the activity trend of the cycled samples ( Supplementary Fig. 10b). To directly visualize the reconstruction degree in CV and compare the activities of the pristine and cycled samples within the same context, the current densities of the cycled samples are still normalized by the surface area of the precatalysts. By plotting the experimental activity enhancement (the change of current density @ 1.62 V vs. RHE, Δ I ) and the calculated metal-oxygen covalency polarity of spinel (the metal-oxygen covalency difference between MO 4 and MO 6 , Δ C ) in the same pattern, a reasonably similar trend is present (Fig. 3i). The fitting between these two parameters shows a R square value of 0.99 (Fig. 3i insert), which strongly demonstrates that the metal-oxygen covalency polarity is reasonable to explain, and therefore can be used to predict, the surface reconstruction ability of spinel oxides. In other words, a structurereconstruction relationship of spinel oxides can be outlined, i.e., higher metal-oxygen covalency polarity between MO 4 and MO 6 can result in stronger surface reconstruction, and vice versa.
Guidelines for the design of highly active oxyhydroxides from spinel pre-catalysts Although the profile of surface reconstruction on spinel oxides has been outlined, it remains ambiguous how to maximize the function of surface reconstruction towards optimal reaction performance. This necessitates the understanding of how surface reconstruction contributes to the measured OER activity and what the key factors dominating the activity of oxyhydroxides are. Generally, the activity enhancement of electrochemical reactions can be realized via two aspects, i.e., increasing the number of active sites or increasing the intrinsic activity of active sites 3 . To figure out how surface reconstruction contributes to the OER activity, electrochemical analysis of cycled Li x Co 3−x O 4 is performed. As shown in Fig. 4a, the measured OER activity significantly rises with the increase of amorphization degree, indicating the reaction performance is positively correlated to the degree of surface reconstruction. The intrinsic activities of the generated cobalt oxyhydroxides are compared in Fig. 4b. As shown, the Tafel slopes of the generated cobalt oxyhydroxides are quite close to each other, in which the maximum difference is less than 10 mv dec −1 . This suggests all the cobalt oxyhydroxides show similar reaction mechanism and share the same rate-determining step 25 . The activity per cobalt site is evaluated by calculating the turnover frequency (TOF). As given in Fig. 4b, the TOF values of the five cobalt oxyhydroxides are within the range of 0.01 s −1 to 0.02 s −1 . Considering the observation that TOF can range across many orders of magnitude, these values indicate the intrinsic activity of each cobalt site in the five oxyhydroxides is almost identical to each other 25 . Therefore, we can rationally deduce that, for pre-catalysts with the same metal cation, an increased reconstruction ability does not alter the intrinsic catalytic activity of the active sites in the in situ generated amorphous oxyhydroxide, but contributes to the enhanced activity by increasing the number of available active sites. The demonstrated structure-reconstruction relationship also allows us to design spinel pre-catalysts with the same reconstruction ability, based on which we can in situ generate oxyhydroxides with the same amorphization degree to investigate the influence of cation type and ratio on the OER activity of oxyhydroxides. Specifically, the reconstruction ability of spinel Zn[Al 1.5 M 0.5 ]O 4 (M = FeCo, FeNi, FeMn, CoNi, CoMn, and NiMn) were screened to search for spinel precatalysts with the same reconstruction ability. The cations of iron, cobalt, nickel, and manganese were studied because they are earthabundant and therefore exhibit cost advantage for potential commercial usage. The screening results are summarized in Fig. 4c. It can be seen that, in zinc aluminum spinel substrate, the introduction of octahedral nickel cation and manganese cation generally results in higher and lower reconstruction ability, respectively. The Δ C value of 3.5 eV is chosen as a reference line to select spinel pre-catalysts with the same reconstruction ability. The Δ C value of 3.5 eV is chosen because it guarantees a high-level reconstruction ability of spinel pre-catalyst and can incorporate sufficient screened samples with the same reconstruction ability to design control experiment for the catalytic investigations on the in situ generated oxyhydroxides. Specifically, spinel  Figure 12). Electrochemical measurements were thereafter performed to in situ generate the corresponding amorphous oxyhydroxides. To guarantee a complete reconstruction, the measurements were conducted for at least 500 cycles until a stable electrochemical curve appears. The TEM images of the samples after reaction are shown in Supplementary Figure 13, which evidence the generation of amorphous layers and nearly the same amorphization degree. Besides, the point scan results (Supplementary Figure 14) confirmed the metal-metal ratio in the amorphous layers is about the same as that in the spinel pre-catalysts.
The electrochemical activity of the generated oxyhydroxides were analyzed and shown in Fig. 4d. The results demonstrate the activity of oxyhydroxides can be largely influenced by the cation type and cation ratio. The comparison between Ni 0.25 Fe 0.75 OOH and Ni 0.25 Co 0.75 OOH suggests that the reactivity of iron cation in oxyhydroxide is higher than that of cobalt cation, consistent with the conclusion reported in other reports 11,25 . The comparison between iron-cobalt oxyhydroxides with different Fe/Co ratio, however, shows a trend that the activity increases with the increase of cobalt ratio in oxyhydroxides. This is a phenomenon that inevitably occurs during OER, in which small amount of iron (normally less than 40%) can significantly raise the activity of cobalt-based oxyhydroxides 12,[25][26][27][28] . There are two possible reasons. One is that the CoO 6 units serve as a conductive and chemically stable host for iron species, which are intrinsically very active but are electrically insulating 25,28,29 . The other is that the formation of high-valent iron species (normally refers to Fe 4+ ), serving either as direct catalytical centers or as redox activators for cobalt cations, is the key for the superior activity 27,28 . To identify the dominant factor for the iron-cobalt oxyhydroxides generated from spinel pre-catalysts, we further investigated the charge transport ability and cation oxidation state.
As elucidated elsewhere, the oxygen hole density is an effective index to estimate the charge transport ability within the transition metal oxides (TMOs), where a higher hole density suggests a higher charge transport ability [30][31][32] . The calculated oxygen hole density is presented in Fig. 4e. For FeOOH, only traceable hole density (0.002 μ B ) is observed, indicating an insulating nature. When 25% cobalt cations are introduced (Fe 0.75 Co 0.25 OOH), the oxygen hole density dramatically increases and reaches 0.07 μ B , suggesting a significant increase of charge transport ability. When the cobalt ratio further increases to 50% (Fe 0.5 Co 0.5 OOH), the oxygen hole density undergoes another dramatic increase to reach 0.13 μ B . With further increasement of cobalt ratio to 75%, the oxygen hole density of Fe 0.25 Co 0.75 OOH is 0.14 μ B and is close to that of Fe 0.5 Co 0.5 OOH, which implies that the charge transport ability within iron-cobalt oxyhydroxides may not undergo further increasement with additional incorporation of cobalt cations. Noting that the oxygen hole density on CoOOH is 0.129 μ B , the values of 0.13 μ B on Fe 0.5 Co 0.5 OOH and 0.14 μ B on Fe 0.25 Co 0.75 OOH indicate that their charge transport ability has already reaches the level at which it is no longer the bottleneck of the reaction. Therefore, we may conclude that the poor charge transport ability is the main catalytical drawback of iron-cobalt oxyhydroxide when the cobalt ratio is less than 50% (iron ratio higher than 50%).
The calculated cation oxidation states are shown in Fig. 4f. The two dashed lines represent the oxidation state of cobalt and iron in CoOOH and FeOOH, respectively. The trend evidences a gradual increase of iron oxidation state with the increase of cobalt ratio, which implies the presence of high-valent iron cations in low-iron-containing iron-cobalt oxyhydroxides. Considering the charge transport ability is no longer a catalytical drawback for iron-cobalt oxyhydroxides with cobalt ratio more than 50%, the activity enhancement with additional cobalt incorporation may attribute to the increased iron oxidation state that facilitates an in situ generation of high-valent iron cations under OER potential. Therefore, the design of highly active and cost-effective oxyhydroxides from spinel pre-catalysts should consider two aspects: (i) the high reconstruction ability that guarantees the availability of more active sites and (ii) the incorporation of conductive and high-valent iron cations that guarantees high intrinsic activity of the active sites.

Potential of surface reconstruction ability of spinel oxides
As demonstrated above, the structural origin of surface reconstruction on spinel oxides has been identified and a criterion evaluating the surface reconstruction ability of spinel oxides is unveiled (Fig. 5). To design spinel oxide with high reconstruction ability for the generation of more available active sites, strategies should be made to weaken the metal-oxygen covalency of MO 4 and strengthen that of MO 6 . By contrast, if a spinel pre-catalyst is believed to hold remarkable intrinsic activity for OER, approaches are therefore needed to balance the metal-oxygen covalency between MO 4 and MO 6 to guarantee a stable in situ spinel crystal frame for long reaction durability.
Employing the computational database established in our previous work 17 , a computational high-throughput figure can be plotted to evaluate the reconstruction ability of spinel oxides, which is shown in Fig. 6. The x-axis presents a parameter that can be used to evaluate the possibility of existence of spinel oxides, which is demonstrated in our previous study 33 . According to our findings, the value of O p states relative to M O d states of phase-pure spinel oxides is normally below 1 eV 33 . It can be seen that most of the data points are located within the phase-pure region, ensuring the effectiveness of using this database to screen the properties of spinel oxides. As labeled, spinel LiCo 2 O 4 exhibits tremendous potential for in situ surface reconstruction, indicating it can be used as an ideal pre-catalyst with high reconstruction ability to generate cobalt-oxyhydroxide-host for active cation (such as iron) substitution. Besides, it's worth to note that the Δ C value of spinel LiMn 2 O 4 is only −1.373 eV, which indicates strong resistivity towards reconstruction. This explains why lithium cannot leach out from LiMn 2 O 4 lattice and therefore the spinel structure remains almost unchanged during OER. The overall Δ C range of spinel oxides is from −6 eV to 12 eV, which demonstrates a huge structural flexibility to manipulate the reconstruction degree. Therefore, we can conclude that spinel oxides are a big family of transition metal oxides with high flexibility to serve as pre-catalysts for electrochemical water oxidation.
In summary, combining computational and experimental analysis, we have pinpointed and demonstrated the structure-reconstruction relationship of surface reconstruction on spinel oxides. The metaloxygen covalency polarity of the M T -O-M O backbone is proved to be the structural origin that triggers surface reconstruction, in which higher polarity corresponds to greater degree of reconstruction. With the established structure-reconstruction relationship, the reconstruction ability of spinel pre-catalysts can be precisely predicted and subtly manipulated. Besides, by conducting mechanistic study on a series of oxyhydroxides generated from spinel pre-catalysts, we have showcased how surface reconstruction contributes to OER activity and how activity of the generated oxyhydroxides evolves with the change of cation type and cation ratio. Our work unveils the key fundamentals of surface reconstruction on spinel oxides and provides firm theoretical foundations for the spinel pre-catalyst design. The research scenario 17 The x-axis is a parameter that can be used to evaluate the structure formation tendency of spinel oxides, which was proposed in ref. 33. for studying the surface reconstruction may also be extended to other types of TMOs.

Density functional theory (DFT) calculations
The DFT calculations were performed by the Vienna Ab-initio Simulation Package 34,35 (VASP), employing the Projected Augmented Wave 36 (PAW) method. The exchange and correlation effects were described by using the Perdew-Burke-Ernzerhof 37 (PBE) functional. The cutoff energy was set to be 450 eV in all the cases. The GGA + U calculations were performed using the model proposed by Dudarev et al 38 . The effective Hubbard U values were set to be 4.5 eV, 4.7 eV, 3.3 eV, 4 eV, 6.4 eV, and 4 eV for tetrahedral Co, tetrahedral Zn, octahedral Co, octahedral Fe, octahedral Ni, and octahedral Mn, respectively. The Monkhorst-Pack grids were set to be 6 × 6 × 6 and 9 × 9 × 9 for performing the structural optimization and density of states calculations, respectively. The conventional standard cell of spinel oxide is employed in this work, in which 8 tetrahedral cations, 16 octahedral cations and 32 oxygen anions are contained. During the model construction, all the lithium cations in the lithium-containing spinels are placed in the tetrahedral sites, as both theoretical 39  The vacancy formation energy was calculated as the energy difference between the vacant spinel together with an oxygen atom and the intact spinel structure, which was defined as where E vacantÀbulk , E oxygen , and E intactÀbulk represent the electronic energy of the spinel with one bulk oxygen vacancy, the chemical potential of one oxygen atom referenced to gaseous oxygen molecule, and the electronic energy of the intact spinel, respectively 42 . Under this definition, a lower value of bulk oxygen formation energy suggests an easier case for generating an oxygen vacancy. The surface transformation energy from spinel to layered oxide was calculated as the energy difference (per surface area) between the reconstructed species and the intact spinel surface, which was defined as E spinel to CoOO = E spinelÀslab À n 1 E CoOO À n 2 Co À n 3 Li where E spinelÀslab , E CoOO , n 1 , n 2 , n 3 , and S represent the electronic energy of the spinel slab model, the layered CoOO species, the number of reconstructed CoOO species, the number of cobalt cations in the tetrahedral sites of spinel slab, the number of lithium cations in the tetrahedral sites, and the surface area of the spinel slab, respectively. Under this definition, a lower value of surface transformation energy suggests an easier transformation from spinel structure to layered structure.

Synthesis of the spinel pre-catalysts
The spinel pre-catalysts were synthesized by a modified conventional sol-gel combustion method. The gel was then decomposed in air at 170°C (with a heating rate of 10°C min −1 to 170°C) for 12 hours to remove the remaining water and ground thoroughly. After applying a further heat treatment in air at 400°C for six hours, phase pure spinel oxide powders were obtained.

Procedures of electrochemical analysis
The OER analysis was operated in a three-electrode cell with a working electrode of glassy carbon flake (effective electrode area: 0.196 cm 2 ), a counter electrode of platinum foil, and a Hg/HgO reference electrode filled with 1 M KOH solution. The 1 M KOH solution was prepared by using electronic grade KOH (purity 99.999%, Aladdin), which excludes the Fe contamination issue 43 . The catalyst electrode was fabricated by the recipe drop-castes method. The catalysts were mixed with acetylene black (AB) at a mass ration of 5:1, then dispersed in isopropanol/water (v/v = 1:4) solvent followed by the addition of Na + -exchanged Nafion as the binder. The mixtures were ultrasonicated for 30 min to reach homogeneous ink. The concentration of oxides in ink is 5 mg/ml, and AB is 1 mg/ml. Before drop-casting, the glassy carbon electrodes were polished to a mirror with α-Al 2 O 3 (50 nm) and washed by IPA and water to clean up completely. Finally, the as-prepared ink (10 ul) was dropped onto glassy carbon electrodes to reach a loading mass of 50 μg ox cm −2 and the electrodes were dried overnight at room temperature.
The pH value of the electrolyte was obtained from pH meter measurement by using Orion StarTM A121 portable pH meter (Thermo Fisher Scientific), from which pH values of 14.00 ± 0.05 were measured in this work. The internal resistance was determined right after the CV measurements by electrochemical impedance spectroscopy (EIS). The EIS measurement was conducted by collecting the AC impedance spectra from 200 kHz to 1 Hz with a voltage perturbation of 10 mV versus a constant potential. The internal resistance was the real part of the impedance where the imaginary part of the impedance was zero in the Nyquist plot. The measured internal resistance is 5.00 ± 0.05 Ω in this work. The OER polarization curves were then corrected to obtain the iR-free potential (E -iR) of the working electrode.

Turnover frequency (TOF) calculations
In this work, the TOFs were calculated via the redox method, which is commonly employed in the system of oxyhydroxides. Specifically, the TOF calculation is based on the following equation: where I ox , A ox , E, e, and N represent the BET normalized current density at the overpotential of 300 mV, the oxide surface area measured by BET, the measured current on the electrode, the charge of a single electron, and the number of cobalt atoms in the electrode. The number of cobalt atoms (N) is estimated by the integration of Co redox peak at ≈1.15 V vs. RHE (assuming a one-electron process of Co 2+ /Co 3+ ).

Derivation of surface reconstruction degree (D r )
The surface reconstruction degree was evaluated on the basis that the loaded catalysts can be considered as a sphere model. With the measurement of Brunauer-Emmett-Teller (BET), the total surface area of the loaded pre-catalyst can be obtained by m loading ÁBET. So the following equation can be obtained m loading ÁBET = 4πðR loading Þ 2 : Based on above equation, the radius of the loaded pre-catalyst (R load ) can be obtained as The volume of the loaded pre-catalyst can therefore be evaluated as The reconstructed part can be evaluated by the solvated lithium cations in the electrolyte. By conducting inductively coupled plasma (ICP) analysis, the content of the solvated lithium cations (c Li , with the unit of g) can be obtained. The mole of lithium cations can then be obtained as c Li 6:94 , where 6.94 (g mol −1 ) is the molar mass of lithium. The number of lithium cations in the solution can be evaluated as c Li where V mol is the volume of one spinel molecule obtained from DFT calculations. Based on above derivations, the reconstruction degree (D r ) can be evaluated as ð10Þ

Materials characterization
The X-ray diffraction (XRD) of spinel oxides were carried out on Bruker D8 diffractometer at a scanning rate of 2°min −1 , under Cu-Kα radiation (λ = 1.5418 Å). The X-ray photoelectron spectroscopy (XPS) measurements were performed using PHI-5400 equipment with Al Kα beam source (250 W) and a position-sensitive detector (PSD) was used to determine the surface composition of the materials. Fourier transform infrared spectroscopy (FTIR) data were recorded using FTIR Perki-nElmer Spectrum. The inductively coupled plasma-optical emission spectrometry (ICP-OES) characterizations were performed using Agilent 720 with a detection limit of 0.02 ppm. The sample powders were weighed (accuracy 0.1 mg) and put in a dry beaker. Then 4 mL of aqua regia solution was added into the beaker and heated at 80°C until the powders were fully dissolved. The solution was allowed to cool naturally to room temperature in fume hood before being transferred into a 25 mL volumetric flask. Deionized (DI) water was added to dilute the solution to 25 mL and was shaked thoroughly for ICP-OES testing. The X-ray absorption near edge spectroscopy (XANES) and extended X-ray absorption fine structure (EXAFS) experiments were conducted in transmission mode at beamline XAFCA of the Singapore Synchrotron Light Source 44 . The data reduction and data analysis were performed with the Athena, Artemis and IFEFFIT software packages. The EXAFS data were extracted from the measured absorption spectra using standard procedures. A cubic-spline-fit procedure was used to determine the post-edge background with a subsequent subtraction. Normalization was performed by dividing the data by the height of the absorption edge at 50 eV. The results of Brunauer-Emmett-Teller (BET) measurements are summarized in Supplementary Figure 2 and Supplementary Table 3.

Data availability
The data supporting the findings of this study are available within the article and its Supplementary Information. Additional data are available from the corresponding author upon reasonable request.